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Abstract 

We analyze the formation of polarization-squeezed light in a medium with 
electronic Kerr nonlinearity. Quantum Stokes parameters are considered and 
the spectra of their quantum fluctuations are investigated. It is established 
that the frequency at which the suppression of quantum fluctuations is the 
greatest can be controlled by adjusting the linear phase difference between 
pulses. We shown that by varying the intensity or the nonlinear phase shift per 
photon for one pulse, one can effectively control the suppression of quantum 
fluctuations of the quantum Stokes parameters. 
PACS: 42.50.Dv, 42.50.Lc. 

Keywords: ultrashort light pulse, self- and cross-phase modulation, electronic Kerr non- 
linearity, quantum Stokes parameters. 

I. INTRODUCTION 

Knowledge of the quantum properties of the polarization of light is essential when treating 
issues concerning Einstein-Podolski- Rosen (EPR) paradox [1] or Bell's inequality [2]. The 
quantum analysis of the light polarization is generally based on Stokes parameters associated 
with Hermitian Stokes operators [3,4]. We call a quantum state to be a polarization-squeezed 
(PS) state if it has the level of quantum fluctuations of one of the quantum Stokes parameters 
smaller than the level corresponding to the coherent state [5]. Such a state can be generated 
by using the optical parametric amplification (x 2 ) [6,7]. Recently, effective methods of PS 
state formation using Kerr-like media (x 3 ) and optical solitons in fibers have been proposed 
[5,8-11]. In the Kerr medium the self-phase modulation (SPM) occurs, leading to quadrature 
squeezing with preservation of the photon statistics [12]. 

It was noted in [13] that a quantum treatment of SPM of ultrashort light pulses (USPs) 
must account for the additional noise related to the nonlinear absorbtion. The study of the 
nonlinear propagation in Raman active media [14] accounts for the quantum and thermal 
noises as a fluctuating addition to the relaxation nonlinearity in the interaction Hamiltonian. 
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However, if we deal with USPs propagation, for instance through fused-silica fibers, the 
electronic motion on ~ 1 fs time scale contributes with about 80% to the Kerr effect while 
Raman oscillators give only 20% contribution [15]. An attempt to develop the quantum 
theory of pulse SPM has been undertaken in [16], where the electronic Kerr nonlinearity is 
modelled as a Raman-like one. The electronic nonlinearity was considered correctly with 
the interaction Hamiltonian [17] and the momentum operator [18] in the normally ordered 
form. 

If two-mode radiation with orthogonal polarization and/or different frequencies prop- 
agates through a Kerr medium, then the cross-phase modulation (XPM) and parametric 
interaction can also occur. The approach in [17,18] has been extended for the combined 
SPM-XPM of USPs in [19] , by considering the spectra of fluctuations of quadratures. No- 
tice that, for interacting solitons in a Kerr medium with anomalous dispersion (k < 0), the 
XPM induces the transient photon-number correlations while SPM answers for the intrapulse 
ones [20]. 

The experimental realization of the polarization squeezing implies the spatial overlapping 
of an orthogonally polarized strong coherent beam with squeezed vacuum on a 50 : 50 
beamsplitter [21] or the interference of two independent quadrature-squeezed USPs produced 
in a fiber Sagnac interferometer [8]. The last series of experiments [10,11] involve both the 
temporal and spatial overlaps of two orthogonally polarized quadrature- squeezed USPs. 
From the theoretical point of view, the quantum fluctuations of Stokes parameters in the 
case of the nonlinear propagation in an isotropic Kerr medium were introduced by the 
study of the light polarization [4]. The time-independent quantum treatment of two-mode 
interaction [5] predicted the possibility of forming the PS state in anisotropic Kerr media. 

In this letter we report on the theoretical computation of the spectra of fluctuations of 
the Stokes parameters in the case of two orthogonally polarized propagating USPs in an 
anisotropic electronic Kerr medium. We employ simultaneously the SPM for quadrature 
squeezing and the XPM for the PS state formation and its control. The estimations are 
done in the frame of the quantum theory for SPM-XPM of USPs developed in [19]. The 
response time of the electronic Kerr nonlinearity is accounted and the dispersion of linear 
properties is described in the first approximation of the dispersion theory. We begin with a 
brief description of the model introduced in [19] and based on the momentum operators for 
the pulse fields. For consistency, we present some elements of the algebra of time-dependent 
Bose operators. By using such algebra we get the average values, the correlation functions, 
and the spectra of quantum fluctuations of Stokes parameters. Our concluding remarks close 
the paper. 

II. QUANTUM THEORY OF COMBINED SPM-XPM EFFECT IN 
ELECTRONIC KERR MEDIUM 

The quantum theory developed in [19] is based on the following momentum operators: 
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where j is the pulse number, (3j and (3 are nonlinear coefficients = /icuo^o^O) / SnoV 
[19,22]) (j = 1, 2), N is the normal ordering operator, H{t) is the nonlinear response function 
[H(t) ^ at t > and H(t) = at t < 0], rij(M) = ^"(*, z) is the "photon 

number density" operator in the cross-section z of the medium, and Aj(t,z) and A^(t,z) 
are the photon anihilation and creation Bose operators with the commutation relations 
[Aj(ti, z), A~l(t 2 , z)} = Sj k S(t 2 — tx). The thermal noise is neglected and the expressions 

(1) and (2) are averaged over the thermal fluctuations. We assume that the Kerr effect 
is produced by the electronic motion with characteristic time r r < 1 fs. In our approach 
the pulse duration r p is much greater than the relaxation time r r , and the Kerr medium is 
lossless and dispersionless; that is, the pulses frequencies are off resonance. 

In general, the induced third-order polarization P^ of the medium at the frequency cu 
is: T (w ) = £ j [xii22 Ei (u)Ej (u) E* (u) + X1212 Ej {u)Ei(u)E* (u) + xfLEMEME*^)] 
[5,23]. Since the first two terms in P^ involve operator products accounted in (1) and 

(2) , the last term introduces products of the form Af A\ + A\ A\ which are neglected 
since they are connected with the parametrical interaction. The conditions (also assumed 
here) under which such terms can be neglected were discussed in [5], where it was noted 
that the parametric interaction can be neglected if Az ^> 1 and the phase mismatch A = 
2(k 2 - h) is large, namely |A| > min^A^O)!, 0|A 2 (O)|}. Here = 8mu 2 ( X u22 + X1212 + 
Xi22i)/c 2 (^i + k 2 ), and kj = (uj/cjrij. Note that for circularly polarized modes the last term 
in p( 3 ) is absent, which is valid in isotropic Kerr media where modes obey only SPM and 
XPM. For elliptically polarized fields the last term in induces the circular birefringence 
An c = — (27r/n )2xi22i(l-^+( c,J )| 2 — \E-{uj)\ 2 ) [23]. However, for the successful experimental 
overlapping of the two quadrature-squeezed USPs a birefringence compensator is required 
[10,11]. In general, An c can be neglected when USPs propagate along some symmetry 
directions [23]. 

Since in our case the Kerr effect is of an electronic origin, in the absence of one- and 
two-photon and Raman resonances, H(t) can be chosen as H(t) = (l/r r ) exp (— t/r r ) at 
t > [22]. Then with (1) and (2) the space-evolution equation for Ai(t,z) in the moving 
coordinate system (z — z , t — t — z/u, where t is the running time and u is the group 
velocity), 

- fa® + 2 (t)]Mt, z) = 0, (3) 

has the solution given by 

A 1 (t,z) = eP 1 ® +d *®A 0tl {t). (4) 

Here 6j(t) = i7j9["oj(*)]> Oj(t) = ^[^o,j(^)], 7j = PjZ, 7 = Pz, n ,j(t) = i^(t)i O j(0 
is the "photon number density" operator at the entrance into the medium z — 0, and 
q[no,j(t)] = J^oo h(ti)n j(t — ti)dti, where h(t) = H(\t\). A 2 (t, z) can be derived by changing 
the indexes 1 <-> 2 in (4). In comparison with the so-called nonlinear Schrodinger equation, 
used in the quantum theory of optical solitons (see [20] and refs. therein), in (3) the pulse 
dispersion spreading in the medium is absent. This approach corresponds to the first ap- 
proximation of the dispersion theory [22]. Note that the structure of q[h j(t)} is like that 
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of the linear response in the quantum description of a USP spreading in the second-order 
approximation. 

The statistical features of pulses at the output of the medium can be evaluated by using 
the algebra of time-dependent Bose operators [17-19]. In such algebra we have 

i 0j (ii)e^' (t2) = e^+^-^io^ti), i„j(*i)e d ' (t2) = e^+^-^iojfa), (5) 

where T>j(t 2 — ti) = ijj h(t 2 — ti), V{t 2 — ti) = i^h(t 2 — ti), 2 r yj is the nonlinear phase 
shift per photon for the j-th pulse, and 7 is the nonlinear coupling coefficient. By using the 
theorem of normal ordering [17,18], one obtains the average values of the Bose operators 
over the initial coherent summary state \ao(t)) = |ao,i(i)) <8> \ao, 2 (t)): 

( £ dj(t)^ — e i<t>j(t)-^(t) ^ ( e ^'(*)) — e *^(*)-Aj(*) ^ (5) 

^ e Oj{ti)+Oj{t2)^ = e i[^(ti)+^(t2)]-^(ti)+w(*2)]-/C J -(ti,t 2 ) ) (7) 
( e c\(ti)+G%(t 2 )^ = e j[^(ti)+^(t 2 )]-[w(ti)+w(t2)]-^(ti,t2)_ (g) 

The parameters <f>j(t) = 2^jn j(t), fij(t) = T^o,j(*)/2 are connected with SPM of the 
j-th pulse and (j>j(t) = 2^n j(t), p,j(t) = 7 2 n j(t)/2 with the pulses' XPM. Here <j)j(t) 
(4>j(t)) is the nonlinear phase addition caused by SPM (XPM). The eigenvalue a j(t) of 
the annihilation operator A j(t) over the coherent state \aoj(t)) can be written as cxoj(t) = 
I a o,j (t) I e 1 ^ '® , where (fj(t) is the linear phase of the j-th pulse. Then |o;oj(t)| 2 = (no,j(t)) = 
n j(t). For simplicity let n j(t — 0) = n j. We separate the time dependence in | ck ,j (i) ] by 
introducing the pulse's envelope rj(t) so that |aoj(£)| = \ a o,j(^)\ r j{t) with r,-(0) = 1. In (7)- 
(8), K,j{t u t 2 ) = /i 0j r 2 (ti + 7-/2)0(1-) and ^-(ii,i 2 ) = fxojrjfa + r/2)g(r) are the temporal 
correlators, where /i ,j = Ij^oj, /Jqj = 'jnoj, g{j) = (1 + |r|/r r .)/i(r), and r = £ 2 — £i- 



III. QUANTUM STOKES PARAMETERS AND THEIR AVERAGE VALUES, 

POLARIZATION DEGREE 



In classical optics the polarization state is visualized as a Stokes vector on the Poincare 
sphere and is characterized by the four classical Stokes parameters {Sk}k=o,3- Since So 
defines the total intensity of pulse field, {Sk}k=i,3 characterize light polarization and form a 
Cartesian axis system. Each point on the sphere corresponds to a definite polarization state, 
whose variation is characterized by the motion of the point on the sphere. In quantum optics, 
{£fc}fc=o,3 are replaced by the operators {S k } k=0>3 which obey: [S , Si] = 0, [Si, Sj] = 2ie ijk S k . 
We define the quantum Stokes parameters as: S (t,z) = Z)|=i nj(t, z), §i(t,z) = fii{t,z) — 
n 2 (t,z), 

S 2 (t, z) = A+(t, z)A 1 (t, z) + A+(t, z)A 2 (t, z), (9) 
S 3 (t,z)=i [At(t,z)A 1 (t,z)-At(t,z)A 2 (t,z)}, (10) 

where in our case Aj(t,z) is given by (4) and Aj(t,z) is its Hermitian conjugate. By using 
Eqs. (5)-(6) we compute the average values of {Skit, z)}k=o,3 on the coherent summary state 
\a (t)). Finally, we get: (S (t,z)) = T, 2 j=1 n ,j(t), (Si(t,z)) = fi 0A (t) -n 0j2 (t), and 
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(S 2 (t, z)) = 2[n 0il (t)n , 2 (t)] 1/2 e-[ Al W +A2 W] cos [* 2 (f) - ^(t)}, (11) 
z)) = 2Ki(t)n , 2 (t)] 1/2 e- [AlW+A2(t)1 sin [<§ 2 (i) - ^(t)}, (12) 

where Aj(i) = Hj(t)+p,j(t), $j(t) = (f>j(t) — (f>j(t)+(pj(t). By measurements, one obtains a set 
of measurable quantities associated with the operators (9) and (10). However, the presence of 
quantum fluctuations results in an uncertainty for the measured quantities. The fluctuation 
uncertainty of each (Sk{t,z)) can be associated with a particular region of uncertainty on 
the Poincare sphere. In the case of the nonlinear propagation, the ball-region of uncertainty, 
specifically for the coherent state of light, changes to the ellipsoid of uncertainty [5]. 

In classical optics the polarization degree V is the ratio of the intensity of the polarized 
part of the radiation 7 po i to the total intensity J tot , and is connected with the classical 
Stokes parameters by V = TZ/{Sq), where 1Z is the radius of the classical Poincare sphere, 
1Z 2 = J2k=i(Sk) 2 ■ Since for completely polarized light V — 1, for partially polarized light 
< V < 1. The definition of the quantum polarization degree follows the classical one: 
V(t,z) = K(t,z)/{S (t,z)), where TZ 2 (t,z) = ELi(&(*» z)) 2 . With the expressions (11)- 
(12) we have 

V(t, z) = {l- 4n o ,i(0% 2 (0Ki(0 + ^o, 2 (t)r 2 (l - e -2[Ai(*)+A 2 (t)] )} i/2_ (13) 

Since in real-life situations Ai(t) + A 2 (t) <C 1, we have V(t, z) ph 1. Indeed, a recent study 
[3,5] investigated the nonlinear behavior of V(t, z) for the nonclassical states of light and 
revealed that the deviation of V(t, z) from 1 is a pure quantum effect. Note that the results 
(11)-(13) are similar to the ones obtained in [5]. 

IV. CORRELATION FUNCTIONS AND SPECTRA OF QUANTUM STOKES 

PARAMETERS 

Since the Kerr medium is assumed to be lossless and dispersionless, the operators So(t, z) 
and Si(t, z), as well as their dispersions, are conserved. Therefore, we focus our attention to 
the quantum fluctuations of S 2 (t,z) and S 3 (t,z) by defining their correlation functions as 

Rs k (h, h) = (S k (t u z)S k (t 2 , z)) - (Skih, z)) (S k (t 2 , z)) (k — 2, 3). (14) 

The correlation functions (14) can be analytically computed by using Eqs. (5) and (7)-(8). 
Here we simply write down the result for Rs 2 (t, t + r) in the approximation jj, 7 <^ 1: 

R S2 (t,t + r) = 5(t) + h(r)[no tl (t)(fa(t) - noAQM*)] sm2[$ 1 (t) - $ 2 (t)] 

+ g{r){n^{t)[4>l{t) + ~4>^ (15) 

where for simplicity we denote t\ — t, t 2 — t + r. In the absence of SPM [(f>j(t) = 0] and 
XPM [<f>j(t) = 0] of pulses, Eq. (15) gives us the correlation function for the coherent state 
Rg^(t, t + r) — 8(t), as expected. The spectral densities of quantum fluctuations of S k (t, z) 
can be evaluated by using Wiener-Khintchine theorem: Ss k (oo, t) = /f^ Rs k {t, t + r)e luJT dT. 
Allowing for a small change of the envelope during the relaxation time we obtain 

SsMt) = l + 2L(ft)[n 0i i(t)^ 2 (t) -no, 2 (t)0i(t)]sin2[$i(t) - $ 2 (t)] 

+ AL 2 {n){n^{tM{t) + ^]+n,^ (16) 
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where L(Q) = 1/(1 + Q 2 ), and Q = uor r is the reduced frequency. The spectral density 
(16) depends on the relaxation time r r and quasi-statically changes with time. Besides, the 
second term on r.h.s. of Eq. (16) indicates that the quantum fluctuations of S 2 (t, z) can be 
less than those corresponding to the coherent state S c s ° 2 h {Vl,t) = 1. The correlation function 

of Sz(t, z), as well as its spectrum, can be easily obtained by shifting (15) and (16) in phase 
with 7r/2. 

For an arbitrary fl Q = ujqt T) at which the linear phase difference between incoming pulses 

a ,.\ 1 + ( n ,i(t)(j)2(t) - no )2 (0^i(*) 

Aip(t) opt = - arctan 1 



L(fio){no,i(*)[^(0 + 01(01 +n ,2(tM(t) + 4>l(t)}}, 

+Mt)-Mt)-Mt) + Mt), (17) 

is optimized, the expression (16) reaches the minimum value: 

S S2 (n ,t) = 1 + 2L 2 (fi ){n 0il (t)^(t) +%(t)] + n , 2 (t)[cf> 2 1 (t) + $(*)]} 
-2L(n )([n OA (t)<j) 2 (t)-no,2(t)Mt)} 2 

+L 2 (n ){n 0A (tM(t) +%(t)] +n , 2 (tM(t) + 0?(t)]} 2 ) 1/2 . (18) 

To characterize the deviation of Ss 2 (Q,t) from the coherent level, we define the normal- 
ized spectral variance Sg 2 (Q,t) = [Ss 2 (Q,t) — l]/n ,i(t). In the case of the suppression of 
quantum fluctuations —1 < S S2 (Q,t) < 0. Let us investigate the S S2 (Q,t) by choosing the 
linear phase difference between pulses to be optimal at a defined reduced frequency Qq, and 
change the intensity of one pulse (the control pulse) in comparison with the intensity of the 
other one. Such dependence of S s (Q, t) at t — and Q = on the maximum nonlinear phase 
addition o ,i = <f>i(t — 0) in case the linear phase difference is optimal at f2 = and Q — 1, 
is displayed in Figs. 1 and 2, respectively. Thus, at f2 = the increase of the control pulse 
intensity produces a uniform suppression of quantum fluctuation of S 2 (t, z) for any o ,i > 1- 
At Qq — 1 the increase of the control pulse intensity produces the suppression basically in 
the domain o ,i < 1- The normalized spectral variance S s (Q,t) at t — 0, Qq — 0, and 
fixed 0o,i = 2, is presented in Fig. 3. Now the suppression in S 2 (t,z) is maximal at Q ~ 0. 
Summarizing, the choice of the linear phase difference allows us to obtain the spectra with 
the form of interest, and the increase of the control pulse intensity can effectively control 
the suppression of quantum fluctuations of quantum Stokes parameters. 

One can also control the suppression of quantum fluctuations of S 2 (t,z) by increasing 
the nonlinear coefficient 72 in comparison with 7^ This is equivalent to the increase of the 
Kerr electronic nonlinearity for one pulse n 2 ( 2 ) in comparison with the one for another pulse 
n 2(i) (7j = Pjz)- The variance Sg 2 (Q,t) at t = 0, fi = for various relations between 
nonlinear coefficients 71 and 72 is displayed in Fig. 4 in the simplest case of pulses with the 
same intensity. In this case, the increase of 72 in comparison with r y 1 suppresses the quantum 
fluctuations of S 2 (t,z) basically at low frequencies, Q ~ . A similar dependence, but at 
Qo = 1, is shown in Fig. 5. Now, the squeezing takes place essentially at Q ~ 1 (u> ~ l/r r ). 

Note the experimentally obtained squeezing of —2.8 dB in S 2 reported in [10]. The 
calculations in [10] use quantum noise operators and do not account for the finite relaxation 
time of the Kerr nonlinearity. However, the relaxation time, which was accounted here, is 
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FIG. 1. Normalized spectral variance S s (Q,t) as a function of the maximum nonlinear phase 
addition </>q,i at the reduced frequency fi = for initial phase difference A(p(t) which is optimal at 
S7o = 0. Curves are calculated at time t = 0, 71 = 72/4 = 27, and correspond to no,2 = ^o,i/4 (a), 
no,2 = n ,i/2 (b), n ,2 = n ,i (c), n ,2 = 3n ,i (d). 




_Q 4 I 1 I 1 I 1 I 1 I 1 I 1 I 

0.5 1.0 1.5 2.0 2.5 3.0 




FIG. 2. As in Fig. 1 but for ft = 1- 
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FIG. 3. Normalized spectral variance S s (Q,t) at </>o,i = 2 for initial phase difference A<p(t) 
which is optimal at Qq = 0- Curves are calculated at time t = 0, 71 = 72/4 = 27, and correspond 
to n , 2 = rco,i/4 (a), n , 2 = n 0;1 /2 (b), n , 2 = no,i (c), n , 2 = 3n 0j i (d). 




FIG. 4. Normalized spectral variance Sg (fi, t) at <^>o,i = 2 for initial phase difference Ay?(i) 
which is optimal at S7o = 0. Curves are calculated at time t = 0, no,i = «o,2> 7 = 7i/2 and 
correspond to 72 = 271 (a), 72 = 371 (b), 72 = 471 (c), 72 = (d), 72 = 671 (e), 72 = 771 (f). 
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FIG. 5. As in Fig. 4 but for Q = 1 

of a fundamental importance since it determines the level of quantum fluctuations of Stokes 
operators below the level corresponding to the coherent state [see Eq. (16)]. Besides, we 
indicated the optimal strategy for the successful generation of the PS state in the electronic 
Kerr medium. 

V. CONCLUSION 

We investigated the formation of polarization-squeezed light in a nonlinear medium with 
electronic Kerr nonlinearity. The correlation functions and corresponding spectra of quan- 
tum Stokes parameters S 2 and S 3 were considered. We shown that, by adjusting the linear 
phase difference between pulses, the maximum suppression of the quantum fluctuations of 
S 2 or 5*3 can be realized at the spectral component of interest. It is established that the 
increase of the intensity of the control pulse can be employed to suppress the quantum fluc- 
tuations of £>2- We find that the increase of one nonlinear coefficient (72) in comparison with 
another one (71) produces a substantial suppression of quantum fluctuations of S 2 - 
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